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Abstract. The composite stellar populations of galaxies comprise stars of a wide range of metallicities. 
Subsolar metallicities become increasingly important, both in the local universe when going from early 
towards later galaxy types as well as for dwarf galaxies and for all types of galaxies towards higher 
redshifts. 

We present a new generation of chemically consistent evolutionary synthesis models for galaxies of various 
spectral types from E through Sd. The models follow the chemical enrichment of the ISM and take 
into account the increasing initial metallicity of successive stellar generations using recently published 
metallicity dependent stellar evolutionary isochrones, spectra and yields. 

Our first set of closed-box 1-zone models does not include any spatial resolution or dynamics. For a 
Salpeter initial mass function (IMF) the star formation rate (SFR) and its time evolution are shown to 
successfully parameterise spectral galaxy types E, Sd. We show how the stellar metallicity distribution 
in various galaxy types build up with time to yield after '^12 Gyr agreement with stellar metallicity 
distributions observed in our and other local galaxies. 

The models give integrated galaxy spectra over a wide wavelength range (90. OA- 160/im), which for ages 
of 12 Gyr are in good agreement not only with observed broad band colours but also with template 
spectra for the respective galaxy types. 

Using filter functions for Johnson-Cousins U,B,V,Rc,Ic, as well as for HST broad band filters in the 
optical and Bessel & Brett's NIR J,II,K filter system, we calculate the luminosity and colour evolution 
of model galaxies over a Hubble time. 

Including a standard cosmological model (Hg = 65, JIq = 0.1) and the attenuation by intergalactic hy- 
drogen we present evolutionary and cosmological corrections as well as apparent luminosities in various 
filters over the redshift range from z ~ 5 to the present for our galaxy types and compare to earlier models 
using single (=solar) metallicity input physics only. We also resent a first comparison of our cc models 
to HDF data. A more detailed comparison with Hubble Deep Field (HDF) and other deep field data and 
an analysis and interpretation of high redshift galaxies in terms of ages, metallicities , star formation 
histories and, galaxy types will be the subject of a forthcoming paper. 
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1. Introduction 

The number of high and very high redshift galaxies is 
increasing rapidly these days. Thousands of U- and B- 
dropout galaxies with photometric redshift estimates are 
known from the HST Hubble Deep Fields - North and 
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South (WiUiams et al. 1996, 1998, Gardner et al. 2000, 
Ferguson et al. 2000) - and from ground-based deep sur- 
veys like the VLT- Fors Deep Field (FDF) and others. 
Hundreds of these have spectroscopically confirmed red- 
shifts up to z ~ 6 (e.g. Hu et al. 1999). Deep surveys are 
being conducted at all wavelengths from UV through IR 
and far into the sub-mm range. These brilliant data re- 
quire mature galaxy evolution models for adequate inter- 



pretation. Ideally, these models should cover all the ob- 
servational wavelength baseline to allow for a consistent 
interpretation of all the available data, be as comprehen- 
sive and realistic as possible, and extend from very early 
phases at very high redshift towards the global properties 
of nearby galaxies of various spectral types. Moreover, an 
ideal galaxy model should be as simple as possible, in- 
volving the smallest possible number of free parameters. A 
comprehensive galaxy model should describe the evolution 
of as many observable quantities as possible (spectrum, 
luminosities, colours, emission and absorption features for 
the stellar population, the gas content and a large number 
of element abundances for the interstellar medium (ISM)). 
A realistic galaxy evolution model should consistently take 
into account both the age and metallicity distributions of 
the stellar populations that naturally result from any ex- 
tended star formation history (SFH). 

This is what wc attempt with our chemically consis- 
tent spectrophotometric, chemical and cosmological evolu- 
tionary synthesis model. The chemical evolution aspects 
of this model were presented by Lindner et al. (1999), 
based on earlier less complete stellar evolutionary tracks, 
in comparison with and interpretation of the observed red- 
shift evolution of Damped Lya Absorber (DLA) abun- 
dances. Here, we present the spectral and spectrophoto- 
metric aspects of our unified chemical and spectral evolu- 
tionary synthesis code including current stellar evolution- 
ary isochrones. 

Spectrophotometric and cosmological evolutionary 
synthesis models generally applied in current interpreta- 
tions of high redshift galaxy data arc using solar metallic- 
ity input physics only together with specific parametri- 
sations for the star formation (SF) histories of various 
spectral types (e.g. Bruzual & Chariot 1993, Bressan et 
al. 1994, Guiderdoni & Rocca - Volmerange 1987, 1988, 
Fioc & Rocca - Volmerange 1997, Poggianti 1997). Broad 
stellar metallicity distributions extending down to fairly 
low [Fe/H] are observed in local galaxies of various types. 
Wc will show in how far the consideration of a realistic 
metallicity distribution among the stars will affect model 
predictions for high redshift galaxies. 

The first attempts to account for non-solar abundances 
and their impact on the photometric evolution of galaxies 
go back to Arimoto & Yoshii (1986). Einsel et al. (1995) 
used more recent and complete stellar evolutionary tracks 
and colour calibrations for initial stellar metallicities 10~^ 
to 4 • 10~^ to describe in a chemically consistent way the 
photometric evolution of galaxy types E through Sd. In 
Moller et al. (1997) we introduced the concept of chem- 
ical consistency into the spectrophotometric evolution of 
galaxies using Kurucz's (1992) model atmosphere spec- 
tra for a range of stellar metallicities and investigated 
the time evolution of ISM metallicity and luminosity- 
weighted mean stellar metallicities in various wavelength 
bands. For models that well agree with observed template 
spectra (Kennicutt 1992) of various types (E, Sb, Sd) we 
gave decompositions of the total light emitted at wave- 
lengths from U through K in terms of luminosity con- 



tributions from various metallicity subpopulations. This 
clearly showed the considerable widths of the metallicity 
distributions in all 3 galaxy types as well as the impor- 
tance of luminosity contributions from stellar populations 
with subsolar metallicities. Recently Vazquez et al. (2003) 
presented their evolution synthesis code SPECTRAL in 
application to the irregular galaxy NGC 1560 that also 
accounts for the the simultaneous presence of stars of dif- 
ferent metallicities. They use stellar evolutionary input of 
the Geneva group and describe the spectral evolution of 
NGC 1560 on the basis on chemical evolution models by 
Carigi et al. (1999). 

This paper is organised as follows. In Sections 2 and 
3 we present our chemically consistent galaxy evolution 
model and the various pieces of input physics for differ- 
ent metallicities that it uses. Sect. 4 gives a comparison of 
our model colours and spectra with observations of nearby 
galaxies to show that - after a Hubble time - our mod- 
els do reproduce the observed properties of typical local 
galaxy types E ... Sd. In Sect. 5 we present the results 
for our models E, Sb, Sd in terms of spectra at vari- 
ous ages, apparent magnitudes, evolutionary and cosmo- 
logical corrections (including attenuation) in wavelength 
bands UBVRcIc JHK and a series of HST broad band fil- 
ters as a function of redshift. Results are presented for 
the redshift range < z < 4.8 for a cosmological model 
Ho = 65, rJo = 0.1, with a redshift of 5 assumed for galaxy 
formation in both cases. 

We compare our chemically consistent models to mod- 
els using solar metallicity input physics only and present 
a first comparison with HDF galaxies with spectroscopic 
redshifts in Sect. 6. Conclusions arc given in Sect. 7. A 
further set of models which also includes the influence of 
dust in a chemically consistent way will be presented a 
separate paper (see Moller et al. 2001a, b for first results). 
A detailed interpretation of the HDF data with the models 
presented here will be the subject of a forthcoming paper. 

2. Metallicity observations in galaxies 

For quite some time, observational evidence has been ac- 
cumulating for sometimes very broad metallicity distri- 
butions of stars in normal galaxies. Stars in the Milky 
Way disk and halo span a range of at least 4 orders of 
magnitude in metaUicity -4.2 < [Fe/H] < -1-0.3. While 
some years ago, the focus was on super-solar metallic- 
ities e.g. in (the centres of) massive ellipticals, bulges. 
X-ray halos around ellipticals, and the hot intra clus- 
ter medium (ICM), by today, it is clear that the av- 
erage metallicities in all those cases are subsolar. 
The sun, our reference star, stands out in metallicity 
among solar neighbourhood stars. For F, G, K dwarfs the 
[Fe/H] distributions extend from —0.8 to -1-0.4 (Rocha- 
Pinto & Maciel 1998), while B-stars show ([0/H]) = -0.31 
(Kilian-Montenbruck et al. 1994). When averaged over 
1 Re, with Re := effective radius, line strength gradi- 
ents in ellipticals indicate (Z*) (0.5 — 1) • (CaroUo 
& Danziger 1994). For stars in the bulge of our Milky 



Way (Z*) ~ (0.3 - 0.7) • Zq (e.g. McWilliam & Rich 1994, 
Sadler et al. 1996, Ramirez et al. 2000), for the X-ray gas 
halos around elHptical galaxies ASCA observations give 
0.1 < [Fe/H] < 0.7 (e.g. Loewenstein 1999). 

Characteristic HII region abundances (i.e. measured 
at 1 Rc); which give an upper limit to the average gas 
phase abundance, range from Z > Z© for Sa spirals down 
to < |Z0 for Sd galaxies (e.g. Oey & Kennicutt 1993, 
Zaritsky et al. 1994, Ferguson et al. 1998, van Zee ct al. 
1998). Locally, dwarf irregular galaxies have metallicities 
in the range (2 - 30)% Zq (e.g. Richer & McCall 1995). 
The first spectra of Lyman break galaxies at redshifts 
z 3 — 4 have shown that their metallicities, derived from 
stellar wind features, are considerably subsolar, sometimes 
even sub-SMC (Lowcnthal ct al. 1997, Trager ct al. 1997, 
Pettini et al. 2000, Tephtz et al. 2000). Neutral gas in 
damped Lya absorbers observed to z > 4 shows abun- 
dances in the range — 3 < [Zn/H] < (e.g Pettini et al. 
1997, 1999, Lindner et al. 1999). 

3. Chemically consistent galaxy evolution models 

Prom a very principle point of view it is clear that in con- 
trast to single burst stellar populations like star clusters, 
any galaxy with a star formation history extending over 
one to several Gyr, i.e. much longer than the lifetime of 
massive stars, will have a stellar population that is com- 
posite not only in age but also in metallicity - as confirmed 
by the observations cited above. This is what we intend to 
account for in our chemically consistent galaxy evolution 
models. The basic concept of our evolutionary synthesis 
model for galaxies has been described in detail by Fritze 
- V. Alvenslcben & Gerhard (1994), the extended version 
allowing for a chemically consistent (=cc) modelling is 
described by Einscl et al. (1995) for the photometric evo- 
lution and in detail by MoUer et al. (1997), MoUer et al. 
(1999a), Fritze - v. Alvenslcben et al. (1999) for the spec- 
tral and spcctrophotomctric evolution and by Lindner et 
al. (1999) for the chemical evolution. 

In the following we briefly outline the principle of the 
new concept of chemical consistency which we consider an 
important step towards a more realistic galaxy modelling. 
In contrast to single burst single metallicity stellar popu- 
lations like star clusters (Schulz et al. 2002) our chemically 
consistent galaxy evolution model, solving a modified set 
of Tinsley's equations with metallicity dependent stellar 
yields, follows the metal enrichment of the ISM and ac- 
counts for the increasing initial metallicity of successive 
stellar generations, both with respect to the evolution of 
ISM abundances (Lindner et al. 1999) and to the spec- 
tral evolution as presented here. The evolution of each 
star is followed in the HR diagram from birth to its fi- 
nal phases according to stellar evolutionary tracks appro- 
priate for its initial metallicity so that at each timestep 
the distribution of all stars over the HRD is known. The 
evolution of the HRD population is followed with vari- 
ous sets of stellar isochrones from the Padova group for 
five different metallicities from Z = A - 10~^ to 5 • 10~^. 



Stellar subpopulations formed with initial metallicities in 
between two of the 5 discrete metallicities of the Padova 
isochrones are described by two components from the two 
adjacent metallicities. The relative contributions of these 
two components are weighted by the inverse of the loga- 
rithmic differences between the metallicity of the subpop- 
ulation and the metallicities of the adjacent isochrones. 
At any timestep the HRD population is used to synthe- 
sise an integrated galaxy spectrum from a library of stellar 
spectra. This library comprises stellar model atmosphere 
spectra from UV to the IR for all spectral types and lu- 
minosity classes for 5 metallicities (Lejeune et al. 1997, 
1998). The total galaxy spectrum is obtained by sum- 
ming the isochrone spectra, weighted by the star formation 
rate at birth of the stars on the respective isochrone for 
each metallicity and, finally, by co-adding the spectra of 
the various single metallicity subpopulations. Combining 
the spcctrophotomctric time evolution with a cosmological 
model and some assumed redshift of galaxy formation we 
calculate the evolutionary and cosmological corrections as 
well as the evolution of apparent magnitudes from optical 
to NIR for various galaxy types taking into account the 
attenuation of the emitted galaxy light by intervening HI 
(cf. Sect. 3.4). 

3.1. Input physics 

In an attempt to keep the number of free parameters 
as small as possible, our models are calculated as closed 
boxes with instantaneous and perfect mixing of the gas. 
We use isochrones from stellar evolutionary tracks pro- 
vided by Bertelli et al. (1994) (Padova group) in the ver- 
sion from November 1999 that include the thermal pulsing 
AGB phase as described in Sclmlz et al. (2002) . The stellar 
lifetimes for the various metallicities are fully taken into 
account in our description of the ISM enrichment. 

The only stellar library covering the whole range of 
stellar metallicities, spectral types, and luminosity classes 
is the library of model atmospheres by Lejeune et al. 
(1997, 1998), based on the original library of Kurucz 
(1992). Its wide wavelength range from 0.09 to 160000 
nm extending far into the UV allows us to calculate cos- 
mological corrections even in the U band out to very high 
redshift z~5. For stars hotter than 50.000 K, the high- 
est effective temperature of Lejeune et al.'s library, we 
use black body spectra. To follow in detail the chemical 
evolution of the ISM, Lindner et al. (1999) included stel- 
lar yields from Woosley & Weaver (1995) and van den 
Hock & Groenewcgen (1997) for various metallicities. For 
the models in this paper we are not interested in the de- 
tailed chemical evolution. We only use the time evolution 
of the global metallicity Z to know when to switch from 
one isochrone to the next (= more metal rich). 

With this new set of metallicity dependent input 
physics our models are now chemically consistcint both 
with respect to the spectrophotometric and to the chemi- 
cal evolution. 



The various spectral galaxy types of the Hubble 
Sequence of normal galaxies are described by their respec- 
tive appropriate galaxy- averaged star formation histories. 
For spheroidal galaxies (E) models use a Star Formation 
Rate SFR ~ e"*^ with an e- folding time t, = 1 Gyr. 
Following Kcnnicutt (1998) wc assume for the spirals a 
SFR linearly proportional to the gas-to-total mass ratio 
with characteristic timescales for the transformation of 
gas into stars ranging from t, = 4 Gyr for Sa through 
t* > 15 Gyr for Sd spectral types. 

We use a standard Salpeter IMF from lower to upper 
mass limits mi = O.O8M0 to mu ^ 70 Mq, as given by the 
isochrones. The IMF is normalised to a fraction of visi- 
ble mass of FVM = 0.5 to match the mass-to-light (M/L) 
ratios for today's (^12 Gyr old) galaxies. 

A Scalo (1986) or a Kroupa (1993) IMF would produce 
a lower number of massive stars than the Salpeter IMF. 
Lindner et al. (1999) have shown that the metallicity evo- 
lution of spirals is too slow for a Scalo IMF as compared 
to observations. As compared to the Salpeter IMF that we 
use, these other IMFs, whithout adjustment of the SFRs, 
would produce a larger fraction of low metallicity stars 
and, hence, lead to slightly bluer galaxy colours. 

3.2. Filters and calibrations 

We calculate the apparent magnitudes m^, evolutionary 
and cosmological - corrections for the UBV- Johnson & 
RqIc Cousins filter system which is taken from Lamia 
(1982), for the JHK filters given by Bcsscl & Brett (1988), 
and for all broad band HST WFPC2 filters. Magnitudes 
in all filters are calibrated in the VEGAMAG system. On 
the basis of the time evolution of the model spectra we 
provide the evolution in other filter systems, easily calcu- 
lated by directly folding the filter and detector response 
curves with the model spectra. 

3.3. Cosmological model 

In order to compare our models with data of high redshift 
galaxies, we transform the spectrophotometric time evo- 
lution into a redshift evolution with a set of cosmological 
parameters (Ho,Oo,Zf), where Zf is the redshift of galaxy 
formation. 

The age of a galaxy at redshift z = is given by 

to := tgal(z = 0) := tHubble(z = 0) - tHubble(zf) 

After tgai(z = 0) models arc normalised to the observed 
average absolute luminosities Mg of the respective galaxy 
types in Virgo (cf. Sandage et al. 1985a, b). To obtain 
the observable apparent magnitudes mA from the absolute 
magnitudes Ma given by our models, we calculate the evo- 
lutionary (eA) and cosmological corrections (kA) and the 
bolometric distance modulus BDM(Ho,f2o): 

mA(z) = Ma(z = 0, to) + BDM(z) + exiz) + kxiz) 



The cosmological or k - correction kA in any wave- 
length band A describes the eff'ect of the expanding uni- 
verse, which redshifts a galaxy spectrum of local age to to 
some redshift z 

kA(z) :=Ma(z, to)-MA(0, to). 

kA can also be calculated from observed galaxy spectra. 
In this case, the maximum redshift to which this is possible 
depends on how far into the UV the observed spectrum 
extends. Our model galaxy spectra at z = extend from 
90 A through 160 jLtm and, hence, allow for cosmological 
corrections in optical bands up to z ;» 10. 

The difference in absolute luminosities between two 
galaxies at the same redshift but with different ages is 
described by the evolutionary correction 

eA(z) := Ma(z, tgai(z)) - Ma(z, to). 

Evolutionary corrections, of course, cannot be given 
without an evolutionary synthesis model. It will be shown 
in Sect. 5 for which galaxy types at which redshifts evolu- 
tionary corrections become important. 

It is important to stress that both the evolutionary 
and the cosmological corrections do not only depend on 
the cosmological parameters but also on the SFH, i.e. on 
the spectral type of the galaxy. 

In this paper we present cosmological and evolution- 
ary corrections and apparent magnitudes for cosmological 
parameters (Hq, Oq) = (65, 0.1) and the formation of 
galaxies at redshift zf = 5. This cosmology gives a galaxy 
age of about 12 Gyr after a Hubble time and, hence, makes 
sure that the agreement with colours and spectra of nearby 
galaxies at z=0 is given. Any combination of Ho and flo 
which gives a local galaxy age tgai ^ 10 Gyr can be ex- 
cluded because the very red colours of ellipticals can only 
be reached after > 12 Gyr and globular cluster ages also 
are of order 12 to 15 Gyr. In particular, high values of 
Ho > 80, even in combination with low values for CIq, yield 
galaxy ages less than 10 Gyr. 

3.4. Attenuation 

For very distant galaxies the cumulative effect of neu- 
tral hydrogen stochastically distributed along our lines of 
sight, mostly in the form of Lya clouds, significantly atten- 
uates the emitted light at wavelengths shorter than rest- 
frame Lya at 1216 A. From his analysis of a large num- 
ber of lines of sight to distant and very distant quasars, 
Madau (1995) derived a statistically averaged attenuation 
correction of the form att{X, z) which we include into the 
cosmological corrections of our models. The effect of at- 
tenuation becomes visible in U at z > 2, in B at z > 2.5, 
in V at z > 3.5, etc. 

3.5. Discussion of models 

Our models describe the global spectrophotometric evolu- 
tion of field galaxies. In particular, our E model represents 



the classical model for a normal elliptical galaxy with av- 
erage luminosity and metallicity. A hierarchical or major 
merger origin of ellipticals is not investigated here. The 
closed-boxes we assume for our spiral models clearly are a 
poor approximation only motivated by our intend to keep 
this first models as simple as possible an the number of free 
parameters to a minimum. For a hierarchically accreting 
spiral the total SFH of its ensemble of subclumps on rea- 
sonable time averages shoiild not differ mutch from what 
our simplified models assume in order to get the correct 
colours and spectra at to- The same had been to be true 
for the chemical evolution of an ensemble of subclumps as 
compared to our simplified closed-box models by Lindner 
et al. (1991). The colour evolution hence is not significant 
affected by our simplification and, as far back as the mass 
of a realistic accreting galaxy is of order > 30% of its mass 
today the differences in apparent luminosities can be es- 
timated to be smaller than differences between different 
galaxy types. In a previous paper we have shown that our 
models well describe the stellar metallicity as observed 
by means of absorption indices (MoUer et al. 1997) and 
the chemical evolution of nearby and high redshift spiral 
galaxies (Lindner et al. 1999). 

In this paper we restrict ourselves to models without 
dust. The effect of dust absorption is analysed in MoUer 
et al. (1999b) with our old track based models. A set of 
models including dust will be presented by Moller et al. 
(in prep.). 



4. Comparison with nearby galaxies 

4.1. Colours 

Star formation histories for our model galaxies are cho- 
sen such that after a Hubble time or, more precisely, after 
the evolution time of a galaxy from its formation at red- 
shift Zf to the present at z = as given by the cosmolog- 
ical model, the colours of our model galaxies agree with 
those observed for nearby galaxies of the respective type 
by Buta et al. (1995), e.g. (B - V) = 0.92, 0.61, and 0.50 
for E, Sb and Sd spectral types, respectively. As a conse- 
quence, however, our V — K colours are then bluer by up 
to 0.3 mag for Es and by up to 0.7 mag for late-type spi- 
rals as compared to the observations of Aaronson (1978). 
This is partly due to the fact that our models describe the 
integrated colour of the entire stellar population, while 
the observations refer to the inner parts of the galaxies. 
Taking into account typical observed gradients in V — K 
(c.f. Fioc & Rocca - Volmerange 1999) brings our V — K 
colours into better agreement with observations, but there 
still remains a difference of about 0.5 mag for late-type 
spirals. As we will show in Moller et al. {in prep.), the 
inclusion of reasonable amounts of dust will bring mod- 
els into very good agreement with observations over the 
entire wavelength range from UV through K. 
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Fig. 1. Comparison of our model spectra (solid) at an age 
Tgai ~ 12 Gyr (E, Sa, Sb, Sc) and Tgai ~ 4 Gyr (Sd) with 
Kennicutt's template spectra (dotted). 
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Fig. 2. Time evolution of Sa model spectra for ages 1, 3, 
6 and 12 Gyr (bottom to top). Fluxes are in units of erg 



cm ^ s 
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4.2. Spectra 

In Fig. ^ we compare our model spectra with templates 
from Kennicutt's (1992) atlas for various spectral types. 
Note the very good agreement between our models E, Sa, 
Sb, Sc with ages of about 12 Gyr with the observed spectra 
of NGC 4889 (E4), NGC 3471 (Sa), NGC 1832 (SBb), and 
NGC 6643 (Sc). The spectrum of NGC 4449 (Sd) is best 
modelled by our Sd galaxy with an age of about 4 Gyr in 
agreement with results found by Bruzual et al. (1993) 

We like to point out that the spectral differences 
among galaxies of the same type in Kennicutt's library are 
larger than the differences between the template galaxies 
and our models. 

Note that we have not included the gaseous emission- 
lines or continuum-lines in our models yet. In a next step, 
gaseous lines and continuum emission will be included in 
our cc models of actively star forming galaxies on the basis 
of metallicity dependent Lyman continuous fluxes and line 
ratios as already shown for single burst single metallicity 
models by Anders et al. (2003) 

The complete set of galaxy model spectra for the vari- 
ous spectral types and ages between 4 Myr to 15 Gyr are 
given in machine readable tables and can also be found 
on our homepage http://www.uni-sw.gwdg.de/~galev. In 
FigElwe show the time evolution of model spectra, on the 
example of an Sa model at 1, 3, 6, and 12 Gyr. 

4.3. Metallicities 

Our code follows the evolution of both the average ISM 
metallicity and the luminosity-weighted metallicities of 
the stellar population as seen in different bands. For 
a galaxy age of 12 Gyr, the average ISM metallicity 
is (Zism) = Zq for E galaxies, (Zism) = 1-5 • Zq for Sa, 
(Zism) = 0.8 • Zq for Sb, (Zism) = 0.5 • Zq for Sc, and 
(Zism) = 0.25 • Zq for Sd spirals, respectively. These val- 
ues are in good agreement with observations of the charac- 
teristic {— measured at ~ 1 Reff) HII region abundances 
in the respective spiral types (e.g. Oey & Kennicutt 1993, 
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Fig. 3. Relative luminosity contributions to U, V, K bands 
from stellar subpopulations of various metallicity for 12 
Gyr old E (top) and Sb (bottom) model. 
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Fig. 4. Relative luminosity contributions to U, V, K bands 
from stellar subpopulations of various metallicity for a 6 
Gyr old Sb model. 



Zaritsky et al. 1994, Phillips & Edmunds 1996, Ferguson 
et al. 1998, van Zee et al. 1998). (See also Sect. 2). 

Depending on the SFH of the galaxy, the average stel- 
lar metallicity may differ by various degrees from the ISM 
metallicity and may also be different in different wave- 
length regions where stars of various masses, ages, and 
hence metallicities dominate the light (cf. MoUer et al. 
1996). 

In Fig. 13 we show the relative luminosity contributions 
to U, V, and K bands from stars of our 5 different metal- 



licity subpopulations for 12 Gyr old E and Sb models. 
For each wavelength band, the sum of the contributions 
from the different metallicities adds up to 100 %. Note 
the broad stellar metallicity distribution of the E model 
extending from Z=4 • 10~* to Z=0.05 in good agreement 
with observed stellar metallicity distributions in resolved 
nearby ellipticals and bulges (cf. McWilliam & Rich 1994). 

The distribution differs from band to band. E.g., stars 
with low metallicity, e.g. Z = 0.0004, contribute about 3 
times more light to the U- than to the K-band, while 
these relative contributions are reversed for stars of higher 
metallicity. 

For models with different SFHs the stellar metallicity 
distributions are, of course, quite different. For stars in our 
global Sb model e.g., the metallicity distribution does not 
extend beyond Zq, and the differences between the rel- 
ative contributions of a subpopulation of given metallicity 
to different bands are smaller than in the E-model. For 
Sd models, the stellar metallicity distribution is sharply 
peaked at Z < ^ Z© with small differences only between 
the different wavelength bands and also to the ISM metal- 
licity. This is readily understood as a consequence of the 
long star formation timescale in Sd galaxies (cf. MoUer et 
al. 1997 for details). 

Similar to Fig.|2|3, Fig.^shows the relative luminosity 
contribution of stellar subpopulations of different metal- 
licities to the light in U, V, and K emitted by an Sb galaxy, 
now at a younger age of 6 Gyr only. By that age - corre- 
sponding to a redshift z ~ 0.5 in our cosmology -, no stars 
of solar metallicity were present in this type of galaxy. 

While Figs. O and 21 give the luminosity contributions 
of different metallicity subpopulations at ages of 12 and 
6 Gyr, the time evolution of the luminosity contributions 
of stellar subpopulations of different metallicities to the 
V-band is shown in Fig. for E and Sb galaxy models. 

It is seen that the broad stellar metallicity distribu- 
tion of the E model is already established at very young 
ages despite the present-day small stellar age distribution. 
Note in particular that despite its SFR declining rapidly 
on the short timescale of t, = 1 Gyr, our E model - due to 
its realistic stellar metallicity distribution - will differ sig- 
nificantly in its spectrophotometric evolution from that of 
any single metallicity single burst model often used in the 
literature for the interpretation of E galaxy observations 
(see also Vazdekis et al. 1996, 1997). 

Over more than the last 50% of its lifetime, the V- 
band light of the elliptical model is coming from stel- 
lar subpopulations of 4 different metallicities (Z=0.0004 
- 0.02) at roughly comparable rates. The V-light of the 
Sb-model is seen to have been dominated by stars with 
half-solar metallicity during the second half of its lifetime. 
Only very recently, i.e. at ages > 10 Gyr, solar metallic- 
ity stars gained importance while at early stages < 6 Gyr 
only stars with Z < 1/4 Zq where present. 

Fig. m very clearly shows that at earlier evolutionary 
times, as observed in galaxies at high redshift, stars of 
lower and lower metallicities become dominant in the spec- 
tra of all spiral galaxies, though not for the E-model. And 
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Fig. 5. Time evolution of the Imninosity contribution to 
the V-band of the different metallicity subpopulations to 
the total luminosity for the E (top) and Sb (bottom) mod- 
els 



this is, of course, not only true in the V-band shown here, 
but over all wavelengths. 

Fig. El shows the time evolution of the luminosity- 
weighted mean stellar metallicities, as defined in MoUer 
et al. (1998), in different passbands for E, Sb, and Sd 
models. These luminosity weighted mean stellar metallic- 
ities in certain passbands are what metallicity dependent 
absorption features in the respective wavelength range of 
the integrated galaxy spectra are expected to measure. As 
already indicated in Fig. O for an age of 12 Gyr there is 
a significant difference in the luminosity weighted mean 
stellar metallicity in various bands in the E-model and 
this is seen in Fig. El to have persisted for all ages > 7 
Gyr. The K-band shows the highest stellar metallicity of 
< Zk >^ 0.014, more than 50% higher than that seen 
in V: < Zv >^ 0.008. In the rest-frame K the time evo- 
lution of the seen stellar metallicity is strongest. While 
it hat rapidly reached a maximum of < Zk >max^ Zq 
at an age around 3 Gyr it tremendously decreases to a 
present < Zk >~ 1/3 Zq since then. The light emitted 
in rest-frame U at a time around 3 Gyr, corresponding to 
a redshift z < 1.5 is shifted to observer-frame I, so we ex- 
pect to see much higher metallicities in high-z ellipticals 
than in local ones. 

Due to the longer timescales of SF the metallicity dif- 
ferences between different passbands get much smaller to- 
wards later galaxy spectral types, as seen in Fig. (HI and 
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Fig. 6. Time evolution of the luminosity-weighted stellar 
metallicities in U, V, and K (restframe) for different spec- 
tral types E, Sb, Sd (top to bottom). 

their time evolution gets much slower and monotonic. 
Note the scale differences among E, Sb, and Sd models. 

5. Evolutionary and Cosmological Corrections 

All evolutionary and cosmological corrections published so 
far were calculated with models using solar metallicity in- 
put physics only (e.g. Bruzual & Chariot 1993, Poggianti 
1997, Fioc et al. 1997). Also still widely in use are the 
k-corrections from Coleman et al. (1980) extracted from 
observed template spectra, although no evolutionary cor- 
rections are available in this case. The wavelength range 
of Colman et al. observation (1400 - 10000 A) hmits the 
redshift range for their k-corrections to z< 1 in U and to 
z< 2 in the R-band. 



To cope with observations of very high redshift galaxies 
up to z ^ 5, models are needed that include the far UV as 
well as evolutionary corrections. 

5.1. Comparison with solar metallicity models 

In the following we will focus on the differences as a func- 
tion of redshift between our chemically consistent mod- 
els (cc models) and models using solar metallicity input 
physics only (Zq models). To isolate the effects of the 
chemically consistent treatment from those due to differ- 
ences in the codes of various authors and/or the particular 
sources of input physics chosen (Padova vs. Geneva stel- 
lar evolutionary tracks, observed stellar spectra vs. model 
atmosphere libraries, etc.) we ran models with our 
code and the same solar metallicity input physics as in our 
cc models. The comparison between cc and Zq models is 
done for cosmological parameters (Ho , JIq,) = (65, 0.1). 
Note that the attenuation is not included in this compar- 
ison because it dominates the the e-|-k-corrections at high 
redshifts and would mash the metallicity effect. 

In order to obtain agreement, after a Hubble time, 
with observed colours of the respective galaxy types (and, 
hence, with our cc model colours) the SFHs have to be 
slightly different in the models that do not contain the 
- on average - bluer and more luminous contributions of 
low metallicity stars. The Zq-E model hence is to be de- 
scribed by an exponentially decreasing SFR with an e- 
folding time of t* ~ 2 Gyr as compared to t* ~ 1 Gyr for 
the cc E model. The Z0-Sb model needs t* ^ 9 Gyr as 
compared to t* ~ 7 Gyr for the cc Sb model. And t* > 15 
Gyr for the Zq and the cc Sd model. 

Figs.[7|and|Hlshow a comparison of e- and k-corrections 
in various bands as a function of redshift for the cc and 
Zq models. The most obvious difference between cc and 
Zq models is seen for the E model at short wavelengths. 
While for redshifts z < 1.2, the cc (e -|- k) corrections in 
U and B are positive, making the apparent magnitudes of 
ellipticals fainter, those of the Zq model get increasingly 
negative from z = through z 2.5 in U and to z ~ 3 in 
B, respectively (cf. Fig. d), making ellipticals apparently 
brighter 

The maximum difference for the E models is as much 
as 3 mag in U and 2 mag in the B-band at z ^ 1. 
Hence, cc models with a SFR decreasing exponen- 
tially with an e-folding time of 1 Gyr predict a 
smaller number of ellipticals to be expected at 
z < 2.5 in magnitude limited samples. At redshifts 
> 2.5 the situation is reversed and the cc models appear 
more luminous. But our simple classic E model may not 
applicable at this redshift range. At longer wavelengths, 
e.g. K, differences are very small by redshift 1 and in- 
crease to 1.5 mag by z = 4. 

For z < 2 the differences in U and B between our Sb 
and Sd cc and the respective Zq models are small. At the 
maximum the cc model has a 0.1 mag higher correction. 
At redshifts higher than 2 the differences of the models 
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Fig. 7. Comparison of the (c + k) corrections in cc and Z© 
models for E (top), Sb (middle) and Sd (bottom) galaxies 
in U, B and K band (Ho, rJo,) = (65, 0.1). 



increases and the Zq models have higher corrections than 
the cc models. The difference grow up to 1 - 1.5 mag in 
the U and B band at high redshift. For the K-band the 
Zq models have always an higher e+k-correction than the 
cc models with a maximum difference of 0.4 - 0.5 mag at 
z ~ 2. 

Hence, on the basis of our cc models we expect 
a slightly smaller number of late-type spirals by 
redshifts z < 2 to show up in magnitude limited 
surveys (U and B) as compared to Z© model pre- 
dictions. At z > 2 we expect more late-type spirals 
as compared to the models 

Fig. |S1 presents the decomposition of the (c + k)- 
corrections shown in Fig. |7| into the evolutionary and cos- 
mological corrections eA(z) and kA(z) for the cc and Z© 
Sb model in the U-, B-,and K-bands. 
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Fig. 8. Comparison of cc and Zq models for Sb in the U- 
, B-, and K-band. Decomposition into evolutionary (top) 
and cosmological (bottom) corrections. 



Let us first look at the e-correction. In the Zq model 
B-band e-corrections become increasingly negative from 
z = to z ~ 2 and remain constant around cb ~ —1.5 mag 
all through z > 4. In the cc model ee decreases faster to 
Cb ~ —3 at z ~ 5. Without the strong compensating ef- 
fect of the k-correction the strongly negative CB-correction 
would make the cc Sb model appear much brighter at 
z > 1.5 than the Zq model. At z > 1.5 the galaxies in our 
cosmology have ages < 2.5 Gyr and stars of low metallic- 
ities Z'- 0.0004 - 0.004 dominate the hght in B. The inte- 
grated spectrum still shows a high flux in B so shortly after 
the major SF epoch. The behaviour of the eu -correction 
is very similar to Cb- In K the evolutionary corrections 
are ck ^ —0.5 through z 4 for both, the cc and Zq 
Sb-models. 

For the k-corrections, the situation is quasi reversed. 
Till z 0.5 the kB-corrcction for the cc Sb model increases 
to kB ~ 1 and then it slightly falls to kB ~ 0.4 at z ~ 3. 
Beyond z = 3, kB increases rapidly to 3.7 at z ^ 5. For the 
Zq model kB follows a trend similar to that of the cc model 
but with an offset of ^ —0.3. The ku-corrections show a 
redshift evolution similar to that of ks. However, the ku 
is smaller than kB at lower redshifts (^ 0.4) and nega- 
tive in the redshift range z ^ 1.3 — 2.6. From z ^ 2.6 — 3.5 
the ku-correction increases from to 3 as the Lyman- 
break gets redshifted into the U-band. Beyond z — 3.5 the 
Lyman-break is shifted out of the U passband and hence 
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Fig. 9. K-correction -including the attenuation- for the 
Sb-model. 



the curve flattens. The difference between the cc and the 
Zq model is again an offset of ~ —0.3. The K-band k- 
correction for the cc model is always negative. It drops 
down to kx ^ —0.9 at redshift z ~ 1.3 and rises to kx ~ 
at z = 5. Until z = 1 the difference between cc and Zq 
models is small (< 0.1). At redshifts ^ 1.5 the cc model 
has an offset in kK of ^ 0.3 — 0.4. 

The k-corrections and the e+k-corrections shown be- 
fore do not include the effect of attenuation to highlight 
the differences between the cc and the Zq models. Fig El 
shows the k-corrections for the Sb model including the at- 
tenuation (cf. Fig.ISI). It is seen that the attenuation is the 
dominant effect beyond z = 3 in U and z = 4 in B. 

5.2. Comparison with earlier work 

We compare both our model c- and k-corrections to those 
of Poggianti (1997). She calculated the e and k correc- 
tion for solar metallicity models. Note that her models 
do not include the attenuation by intergalactic HI and 
- due to a shorter wavelength coverage - do not allow 
for cosmological corrections in U for z > 2 and in B for 
z > 2.5. So we compare with our non attenuated cc mod- 
els at this point. We also calculate the e-correction for the 
cosmology given by Poggianti (1997) for this comparison 
(Ho = 50, flo = 0.5). Figs. [TUl and [TTl show the compari- 
son between our E and Sc cc models with the respective 
models form Poggianti (P) in the U- and B-band. Note 
the good agreement of both models. 

6. Redshift evolution of apparent magnitudes 

The redshift evolution of apparent magnitudes in 
U,B,V,Rc,IcjJjII,K and HST filters is given on our home- 
page. 

In Fig. El we show the redshift evolution of appar- 
ent magnitudes for the E, Sa, and Sd model in in HST 
m450-band (-Johnson B) and (m450 - m606) - (B - V) 
colours and compare to the photometric redshift catalogue 
for the HDF galaxies (Sawicki et al. 1997). For each model 
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Fig. 10. Comparison between our cc models and the 
solar models of Poggianti (1997) (P). The top panel 
shows the k-corrections and the bottom panel the e- 
corrections for E and Sc models in U. Note that we 
compare none-attenuated models in Poggianti's cosmol- 
ogy (iJo,f^o) = (50,0.5). 



Table 1. Evolutionary and cosmological corrections for 
(i?o, ^^o) = (65,0.1). The complete version of the table is 
given on our homepage. 
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the evolution is plotted with (thick lines) and without 
(thin lines) attenuation. 

Before we present a very first and preliminary com- 
parison of cc models with LBG data we recall that 1*^' our 
models refer to integrated properties of galaxies and 2"^^ 
our model galaxy types - being described by SFHs or SF 
timescales - refer to a spectral and not to a morphologi- 
cal classification of galaxies. The question how far back in 




1 1.5 
Readshift z 



Fig. 11. Same like Fig.^J but now in the B-band 



time the locahy observed 1-to-l correspondence between 
spectral and morphological galaxy types remains valid is 
another open question. 

The E model becomes apparently fainter at redshifts 
< z < 1 (m450 ^ 24mag), it then starts getting brighter 
again as it approaches its active star formation phase. 
At z 3 the attenuation comes into play and makes the 
model get fainter rapidly. By z = 5 the difference between 
the attenuated and the non-attenuated models increases 
to ~ 10 mag. 

The Sa and Sd models become fainter at redshifts 
< z < 1, with m450 ~ 24mag for Sa and m450 ~ 27mag 
for Sd at z ~ 1. Thenceforward both spiral models get 
moderately fainter until the attenuation comes into play 
at z ~ 3. 

In comparison to the HDP data from Sawicki et al. 
(1997), the Sa model follows the bright limit familiar well. 
The Sd model gets fainter than the observational limit 
of m450 ^ 29mag beyond z = 2 while the classical initial 
collapse E model is too bright at z > 2 to be comparable 
with the data. The lack of luminous galaxies at redshifts 
lower than 0.5 is the selection effect explicitly intended by 
the choice of the HDF. The lack of galaxies as luminous 
as our classical E model, however is real and not due to 
any selection or bias. But on the basis of the luminosity 




Fig. 12. Redshift evolution of B-band luminosity (m450) 
and from the B-V colour (m450-m606) compared to the 
HDF galaxies photometric redshift catalogue of Sawicki 
et al. (1997). Vertical bars show the la luminosity and 
colour ranges of local galaxies 



we can not determine the role of attenuation. To do this 
observations z > 4.5 are needed. 

More information is provided by the colours of galax- 
ies. For the (m450 — m606) colour evolution it is seen that 
for redshifts z > 2 the E model gets as blue as the spi- 
rals. This is partly due to the strong light contributions 
from low metallicity stars and partly due to the youth of 
its stellar population (age ^ 1.6 Gyr at z ~ 2). In the 
redshift range 0.25 ^ z < 0.75 the E model is very red 
(m450 — m606) ^ 1.7. The colour difference between the 
Sa and Sd models is small, ^ 0.3. The models follow the 
data well, but there are a lot of very blue objects which 
cannot be explained by our models of undisturbed galaxy 
evolution. 

At redshifts z ^ 3 the models with attenuation fit 
the data mutch better than models without attenuation. 
While at z ^ 1 a large number of faint and blue galaxies 
are seen in the HDF, bluer and fainter than our late-type 
spiral models, the bulk of galaxies at 1 ^ z ^ 3 are well 
compatible wit our normal spiral progenitor models. 

A detailed modelling and interpretation of these 
Lyman Break Galaxies (LBGs) is beyond the scope of 



the present paper. It will be the subject of a forthcom- 
ing paper including a much larger number of LBGs both 
with spectroscopic and photometric redshifts (see MoUer 
et al. 2001b for first and preliminary results). This investi- 
gation will use the full multi-colour information available 
for HDF and other deep field galaxies. 

7. Conclusions 

Our evolutionary synthesis model is now chemically con- 
sistent with respect to both the spectrophotometric and 
the chemical evolution. This gives the possibility for a de- 
tailed study of the age and metallicity distributions in the 
composite stellar populations of galaxies. In the frame- 
work of simplified closed-box models where the different 
spectral types of galaxies E, Sa, Sb, Sc, and Sd are de- 
scribed by their respective appropriate star formation his- 
tories we analyse how the presently observed average stel- 
lar abundances and abundance distributions have evolved 
in time. 

While - with somewhat different star formation histo- 
ries ~ models using input physics for one single metallicity 
(e.g. solar) only can also be brought into agreement with 
locally observed galaxy colours and spectra, the evolution 
with redshift is significantly different. 

Before analysing galaxy data out to high redshifts 
we made sure that our chemically consistent models cor- 
rectly describe the integrated properties of galaxy types 
E, Sa-Sd in their spectrophotometric appearance from the 
UV through the NIR, their average ISM abundances, the 
metallicity of their stellar populations, their gas content, 
present-day star formation rates and the chemical abun- 
dances of various elements. 

This large number of observables allows to tightly con- 
strain the only free parameter of our closed-box 1-zone 
models: the time evolution of their star formation rate 
or their star formation history. We show how the stellar 
metallicity distribution in various galaxy types build up 
with time to yield after ^ 12 Gyr agreement with stel- 
lar metallicity distributions observed in our and other 
local galaxies. We also presented the time evolution of 
the luminosity-weighted mean stellar metallicities of dif- 
ferent galaxy types in various bands. It is these luminosity- 
weighted metallicities that are expected to be measured 
by metal-sensitive absorption features observed in the in- 
tegrated light in the respective passbands. 

For the spectral galaxy types E, Sa to Sd, we give the 
spectral energy distributions over a wide wavelength range 
(9.09 - 160000) nm in their time evolution from ages of 1 
Gyr up to 15 Gyr. 

Any desired set of filter response functions can directly 
be used to calculate from these spectra the time and red- 
shift evolution of galaxy luminosities, colours, evolution- 
ary and cosmological corrections and apparent magnitudes 
(for any desired cosmological model). 

In comparison with models using solar metallicity in- 
put physics only, we discuss the effects of the inclusion of 
subsolar metallicity stellar subpopulations. In chemically 



consistent models E galaxies appear significantly fainter 
to redshifts z ~ 2.5 as compared to solar metallicity mod- 
els. A smaller number of ellipticals is thus expected to 
contribute to magnitude limited samples. 

Chemically consistent spiral models, on the other 
hand, appear brighter at z > 2 than they would if only 
solar metallicity input physics were used. Hence, we ex- 
pect a larger number of intermediate and late-type spirals 
from z > 2 to show up in magnitude limited surveys. 

We present a large grid of evolutionary and cosmologi- 
cal corrections as well as apparent magnitudes and colours 
in various filter systems (Johnson, HST, ...) from UV to 
NIR including the effect of attenuation by intergalactic 
HI for galaxy types E and Sa to Sd using cosmological pa- 
rameters (Ho, rio) = (65, 0.1) with a redshift of galaxy 
formation assumed to be Zf = 5. 

Models and results our chemically consistent chemical 
evolution models in terms of a large number of individual 
element abundances in the ISM of various spiral types 
were presented in Lindner et al. (1999) and used for the 
interpretation of Damped Lya Absorbers. 

A very first and preliminary comparison with the red- 
shift evolution of HDF galaxies with photometric redshifts 
from Sawicki et al. (1997) indicates that their luminosi- 
ties, colours and SFRs are well compatible with those of 
our normal spiral models over the redshift interval from 
z ^ 0.5 all through z > 4. 

A detailed and extensive comparison of our chemi- 
cally consistent model results with the full set of observed 
colours and luminosities of high redshift galaxies (e.g. 
Lyman Break Galaxies) will be presented in a forthcoming 
paper. 

In the framework of chemically consistent models it 
is also possible to include the effects of dust in a largely 
consistent way, tying the amount of dust to the evolving 
gas content and metallicity. Stellar track based chemically 
consistent models with dust will be presented in a com- 
panion paper by MoUer el al. {in prep., see MoUer et al. 
2001a, c for first results). 

Appendix A: Available data 

The following data are available at our homepage 
http://www.uni-sw.gwdg. de/'^galev/ccmodels/ e- and k- 
corrections, apparent magnitudes, spectra as function of 
redshift, and spectra as function of evolution time given in 
restframe. These data are available for all Models (E-Sd). 
Please see the README file for detailed information and 
file format. 
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